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Recent  analysis  of  direct  numerical  simulations  of  compressible  homogeneous  shear  flow 
turbulence  has  unraveled  some  of  the  energy  transfer  mechanisms  responsible  for  the  decrease 
of  kinetic  energy  production  when  the  flow  becomes  more  compressible.  In  this  complemen¬ 
tary  study,  we  focus  our  attention  on  the  rate  of  strain  tensor.  A  Helmholtz  decomposition 
of  the  velocity  field  leads  us  to  consider  a  solenoidal  and  an  irrotational  rate  of  strain  tensor. 
Their  eigenvalue  distributions,  eigenvector  orientations,  and  the  relative  alignment  between 
the  eigenvectors  and  the  vorticity  and  pressure  gradient  vectors  are  examined  with  the  use  of 
probability  density  functions.  The  irrotational  rate  of  strain  tensor  is  found  to  have  a  pref- 
ered  structure  in  regions  of  strong  dilatation.  This  structure  depends  on  the  mean  shear,  and 
is  quite  different  from  that  of  the  solenoidal  rate  of  strain  tensor.  Compressibility  strongly 
affects  the  orientation  properties  of  the  pressure  gradient  vector. 


*  Research  supported  by  the  National  Aeronautics  and  Space  Administration  under  contract  No.  NASl- 
18605  while  resident  at  the  Institute  for  Computer  Applications  in  Science  and  Engineering  (ICASE),  NASA 
Langley  Research  Center,  Hampton,  VA  23665. 
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1  Introduction 


In  recent  direct  numerical  simulations  of  compressible  homogeneous  shear  flow  [1,  2],  global 
properties  of  the  flow  were  examined  to  determine  the  effect  of  compressibility  on  the  growth 
of  kinetic  energy  and  of  dissipation.  Conclusions  of  this  study  indicate  that  compressibility 
effects  (associated  with  increased  initial  turbulent  Mach  number)  decrease  the  growth  rate 
of  the  kinetic  energy.  This  reduction  in  the  growth  rate  was  determined  in  [1]  to  be  due 
to  enhanced  compressible  dissipation  and  pressure-dilatation  terms,  and  due  to  reduced 
production  by  the  mean  shear  in  the  energy  equation.  The  present  study  emphasizes  the 
properties  of  the  rate  of  strain  tensor  in  order  to  provide  insight  into  the  structure  of  the 
turbulence  at  the  intermediate  and  small  spatial  scales. 

It  is  now  accepted  that  homogeneous  incompressible  turbulence,  although  stochastic  in 
nature,  contains  structures  which  retain  their  identity  over  extended  time  and/or  space. 
These  are  the  so-called  coherent  structures,  and  much  research  is  devoted  to  isolating  and 
understanding  their  kinematical  and  dynamical  properties.  For  example,  the  proper  orthog¬ 
onal  decomposition  (POD)  [3,  4]  constructs  a  set  of  mutual  orthogonal  eigenfunctions  that 
maximizes  the  kinetic  energy  content  of  each  successive  eigenfunction.  Simulations  using 
very  few  of  these  modes  clearly  indicate  that  they  can  capture  the  essential  features  of  the 
flow.  Alternatives  to  the  POD  have  also  been  proposed.  Vorticity  eduction  methods  [5], 
standard  correlation  techniques  [6],  and  wavelet  packet  techniques  [7]  are  but  three  of  them. 
Yet  another  approach  has  been  advanced  by  Ashurst  [8]  who  proposes  to  describe  turbu¬ 
lence  as  a  statistical  ensemble  of  elemental  vortices.  The  geometry  of  the  vortices  and  the 
rules  that  govern  their  interaction  would  then  be  sufficient  to  reproduce  some  of  the  known 
small-scale  properties  of  the  flow. 

Within  the  last  decade,  detailed  analyses  of  numerically  generated  databases  have  re¬ 
vealed  new  structural  properties  of  homogeneous  turbulence,  hitherto  unsuspected.  Most 
of  the  work  was  motivated  by  the  desire  to  understand  the  origins  of  the  intermittent  dis¬ 
sipative  structures.  While  it  has  been  known  for  a  long  time  that  vortex  stretching  is  a 
dominant  mechanism  for  the  production  of  small  scales  [9,  10,  11],  the  theories  were  either 
quite  heuristic,  or  based  on  global  averages.  The  advent  of  powerful  supercomputers  has 
permitted  detailed  visualizations  of  the  flowfields,  and  has  led  to  the  inescapable  conclusion 
that  the  intermittent  structures  are  regions  of  extremely  high  vorticity  and  may  be  classified 
as  coherent  structures  in  the  sense  that  their  lifetime  exceeds  that  of  typical  turbulent  eddies. 
It  is  hoped  that  these  mechanisms  will  offer  deeper  insight  into  the  energy  cascade  process, 
and  thus  lead  to  good  subgrid  scale  turbulence  models.  Some  attempts  in  this  direction  have 
already  been  made  by  Lundgren  [12]  who  proposes  the  concept  of  a  spiraling  vortex.  With  a 
minimum  number  of  assumptions,  he  provides  the  mathematical  framework  from  which  the 
Kolgomorov  power  law  is  deduced. 

Along  with  more  detailed  information  about  the  intermittent  dissipation  regions,  sev¬ 
eral  apparent  paradoxes  were  uncovered  [13,  14,  15,  16].  For  example,  numerical  evidence 
suggests  that  the  vorticity  vector  aligns  itself  primarily  along  the  eigenvector  of  the  rate 
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of  strain  tensor  which  corresponds  to  the  intermediate  eigenvalue  (which  has  a  tendency 
to  be  positive).  This  is  contrary  to  the  expectation  that  the  vorticity  grows  fastest  in  the 
direction  of  maximum  strain,  and  thus  should  be  aligned  in  the  direction  associated  with 
the  largest  eigenvalue.  Further  surprises  result  from  the  presence  of  tubular,  high  intensity, 
vortex  structures  within  which  two  of  the  eigenvalues  of  S,j  are  positive.  One  would  normally 
expect  flat  vortex  structures.  In  an  attempt  to  explain  some  of  these  results,  Cantwell  [17], 
Majda  [18],  Ashurst,  Kerstein,  Kerr  &  Gibson  [13],  and  Vieillefosse  [19]  show  that  some  of 
these  effects  can  be  predicted  from  the  inviscid  equations  if  the  anisotropic  component  of  the 
pressure  Hessian  p,ij  is  neglected.  In  particular,  the  inviscid  equations  predict  the  observed 
positiveness  of  two  eigenvalues  of  Sij,  as  well  as  the  observed  alignment  of  vorticity.  Unfor¬ 
tunately,  this  model  does  not  predict  tne  preferred  shape  (eigenvalues  in  the  ratio  —4:1:3) 
of  the  principal  ellipsoid  of  Sij  discussed  by  Ashurst  et  al.  [13].  This  leads  to  the  hypothesis 
that  the  effect  is  either  related  to  viscosity,  or  to  some  anisotropic  property  of  the  pressure 
Hessian.  With  a  stochastic  model  for  the  viscous  terms,  Girimaji  and  Pope  [20]  reproduce 
many  of  the  statistical  properties  of  the  flow  invariants.  Data  has  also  been  obtained  on  the 
alignment  properties  of  the  pressure  gradient  with  respect  to  the  eigenvectors  of  Sij.  Several 
researchers  agree  that  Vp  aligns  itself  in  the  direction  of  maximum  compressive  eigenvector 
of  Sij,  although  a  plausible  explanation  has  not  yet  been  advanced  [13,  15]. 

This  intrinsic  turbulent  structure,  so  pervasive  in  incompressible  homogeneous  isotropic 
and  shear  turbulence,  is  affected  by  compressibility.  Erlebacher,  Hussaini,  Kreiss  Sarkar  [21] 
showed  that  for  shockless  isotropic  turbulence,  the  dilatational  component  of  the  velocity 
(v^)  decouples  from  the  solenoidal  component  (v^)  to  lov.'est  orderU  More  precisely, 
satisfies  the  incompressible  Navier-Stokes  equations,  which  in  turn  defines  an  associated  in¬ 
compressible  pressure  p^  The  remainder  fluctuation  pressure,  p^,  together  with  satisfy 
a  set  of  equations,  which  only  couple  to  the  nonlinear  solenoidal  flow  through  terms  of  order 
Mt  where  Mt  is  the  (assumed  small)  rms  turbulent  Mach  number.  This  result  is  obtained 
through  a  multiple  time-scale  analysis  [21]  in  which  the  phenomena  associated  with  (v^,p^) 
vary  on  a  0{Mt)  time  scale,  whereas  the  solenoidal  flow  evolves  on  a  0(1)  time-scale.  In 
homogeneous  shear  flow,  with  moderate  to  high  shear  rate  S,  the  situation  is  more  compli¬ 
cated  since  both  and  v*"  now  vary  on  the  same  0(5“^)  time-scale,  and  so  are  expected 
to  interact  more  strongly  with  each  other  than  in  isotropic  flow. 

Very  little  is  currently  known  about  the  small  scale  structure  in  compressible  homo¬ 
geneous  turbulence.  Blaisdell  [22]  has  done  some  statistical  analysis  of  results  from  high 
resolution  direct  numerical  simulations.  A  significant  result  is  that  the  vorticity  and  the 
dilatation  are  statistically  independent  for  all  time  and  are  therefore  decorrelated.  This  is 
borne  out  in  our  own  simulations. 

The  main  objective  of  this  paper  is  to  initiate  a  systematic  investigation  of  the  small 
spatial  structures  present  in  compressible  homogeneous  turbulence.  This  is  accomplished 
through  an  analysis  of  the  rate  of  strain  tensor  calculated  both  from  the  solenoidal  and 
from  the  irrotational  components  of  the  velocity  field.  Statistical  properties  of  the  tensors’ 

^This  decomposition  is  vvell-ciefiiieci  and  unique  in  homogeneous  flows 
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eigenvalues  and  eigenvectors  are  presented.  We  also  look  into  the  relationships  between  the 
pressure  gradient,  the  vorticity  and  the  eigenvectors  of  the  rate  of  strain  tensors,  and  iiow 
these  are  affected  by  the  mean  shear  rate.  At  low  shear  rates,  we  e.xpect  only  slight  departure 
from  incompressible  statistics.  However,  as  the  shear  rate  is  increased,  both  the  solenoidal 
and  the  irrotational  velocity  components  will  vary  in  a  0(6'“')  time,  determined  by  the  mean 
shear.  This  leads  to  possible  interaction  mechanisms  between  the  two  components. 


2  Equations 


The  homogeneous  shear  flow  data  used  for  the  ensuing  analysis  is  obtained  from  the  direct 
numerical  simulation  (DNS)  of  the  full  compressible  .\avier-Stokes  ecpiations  subject  to  as¬ 
sumptions  of  spatial  homogeneity.  The  mean  flow  is  unidirectional  in  the  .r  direction.  <ind 
has  uniform  shear,  5,  in  the  vertical  direction: 

(1) 

The  mean  temperature  and  density  are  uniform.  The  velocity  components  are  e.\i>ressecl  in 
perturbation  form  which  leads  to  coefficients  in  the  eejuations  which  depend  on  the  normal 
coordinate  y.  A  time-d<‘pendent  transformation  is  applietl  to  the  physical  coordinates  to 
replace  this  y  dependence  by  coefficients  which  depend  on  time.  This  transforiiiation  was 
first  used  by  Rogallo  [23]  in  D.NS,  and  is  now  a  standard  tool  to  transform  s./me  standard 
homogeneous  flows  to  a  form  amenable  to  Fourier  spectral  metliods. 

For  completeness,  we  present  the  Navier-Stokes  equations  in  the  form  in  which  they  are 
numerically  solved; 

dtp  +  [pui),i  -  St{pU2),i  =  0  (2) 


dt{pu,)  +  {pujud,j  =  -p,i  +  rp,j  -  Spu 

+  St[pu2Ui),i  +  6  p  —  Str, 2^1  (3) 

dtp  +  UjP,j  +  ipUj^j  =  Stu-zp^i  +  -,StpU2,i  + 

+  (7  ~  l)''t  i'.j'j  —  '2Sll\i  >  +  S^t"  ( 1) 

p  =  pRT  (5) 

where  the  dissipation  function  <I>  and  stres>  tensor  rp  have  their  usual  definition  in  the  labo¬ 
ratory  frame  of  reference.  Other  variables  are  the  fluctuating  velocity  a,,  the  instantaneous 
density  p,  the  pressure  p,  the  tempeiature  T,  the  universal  gas  constant  /?,  the  ratio  of  spe¬ 
cific  heats,  7,  and  the  thermal  ccuiductivity  k.  Derivatives  are  taken  in  a  frame  of  reference 
that  is  moving  with  the  flow. 

Velocity,  density  and  temperature  are  respectively  non-dimensionalized  by  the  initial  vnis 
velocity,  the  constant  mean  density  p,  and  the  constant  mean  temperature,  T.  The  viscosity 
p,  is  constant,  and  the  Prandtl  number  is  set  to  0.7,  valid  for  air. 
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The  coordinate-independence  of  the  flow  equations  in  the  (moving)  computational  domain 
permits  us  to  express  the  solution  in  terms  of  Fourier  basis  functions.  Thus,  the  flow  is 
periodic  and  contained  in  a  cubic  box  of  length  2it.  The  equations  are  solved  using  spectral 
collocation  with  dealiasing.  The  scheme  is  fully  explicit,  and  the  spatial  operator  is  advanced 
in  time  with  a  low-storage  third-order  Runga-Kutta  method.  Details  of  the  algorithm  and 
of  the  initialization  process  are  given  in  [1]. 

Contrary  to  the  dynamics  of  decaying  compressible  Isotropic  turbulence  [21,  26],  in  which 
the  effects  of  initial  conditions  are  never  quite  eradicated,  we  expect  that,  with  proper- 
scaling,  the  structure  of  the  turbulent  shear  flow  is  intrinsic  to  the  flow,  independent  of  the 
initial  data.  The  shear  flow  simulations  are  parametrized  at  initial  time  by  the  root  mean 
square  velocity  Urma  —  y/uiiii,  and  the  degree  of  initial  compressibility,  x  =  The 

superscript  C  refers  to  the  irrotational  component  of  the  velocity  obtained  from  the  (unique) 
Helmholtz  decomposition  of  Ui.  The  solenoidal  component  is  referred  to  by  a  superscript  I. 
Throughout  the  text,  quantities  derived  from  one  or  the  other  of  these  two  velocity  fields 
are  denoted  by  one  or  the  other  superscript.  Further  parametrizations  of  the  initial  field 
include  the  rms  fluctuations  of  two  of  the  three  thermodynamic  variables.  The  density  and 
the  pressure  fluctuations  are  initialized  as  described  in  Sarkar  et  al.  [2]. 

Intrinsic  parameters  useful  to  characterize  the  turbulent  flow  include  the  Taylor  length 
scale  A  =  ^/Urma/ the  microscale  Reynolds  number,  R\  =  and  the  turbu¬ 

lent  Mach  number  M(  =  ^ 


3  Vorticity  and  Dilatation 


While  vorticity  has  long  been  viewed  as  a  fundamental  quantity  in  incompressible  turbulent 
flows,  the  flow  dilatation  is  perhaps  the  single  most  important  measure  of  compressibility. 
The  properties  of  the  dilatation  field  must  be  well  understood  to  extract  the  relevant  lead¬ 
ing  order  compressibility  effects.  The  dilatation,  0  =  u,,,,  describes  to  leading  order  the 
departure  of  the  flow  from  an  incompressible  state.  The  transport  equation  for  0  is 

^  +  iu.“)  -  2S  S.j  +  (/>-■  (6) 

where  Pij  is  the  combined  pressure  and  viscous  stresses. 


Pij  =  -P^ij  +  2/i  Sij  (7) 

0 

and  S ij  is  the  deviatoric  component  of  Sij,  where  Sij  is  the  rate  of  strain  tensor  (u,  j -t-Uj,i)/2. 
The  time  derivative  Df  Dt  represents  the  material  derivative,  following  a  fluid  element.  The 
variance  of  0  satisfies 
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Using  numerical  data,  Sarkar  et  al.  [1]  and  Blaisdell  [22]  find  that  the  dilatation  has  a  prob¬ 
ability  distribution  function  (pdf)  with  negative  skewness.  Although  it  is  clear  that  there  is 
an  inherent  asymmetry  in  a  viscous  flow  due  to  the  entropy  condition,  a  theoretical  deriva¬ 
tion  that  accounts  for  this  skewness  still  eludes  us,  even  when  the  turbulence  is  isotropic. 
The  numbers  beneath  Eq.  (8)  represent  the  budget  of  0  at  St  =  9.  This  data  and  that 
presented  throughout  the  text  correspond  to  DNS  of  homogeneous  shear  flow  with  initial 
conditions  S  =  20,  Mto  =  0.2,  —  1/150,  Rxo  =  13.5  and  incompressible  initial  conditions 

(Prms)o=0  and  xo  -  0.  Production  of  dilatation  results  primarily  from  the  interaction  of  the 

mean  shear  with  0  5^-.  A  lesser  contribution  is  due  the  triple  interaction  term  0  Sfj  Sfy 
The  vorticity  satisfies  the  transport  equation 

=  UJjSij  —  0(u;,  —  SSi^)  —  S  Ssi  +  —  {0J-l8i2  +  U;2<5,i)  +  ^i3k{p~^  Pkp,p),y  (9) 

In  general,  the  viscous  component  can  be  neglected,  and  in  the  absence  of  shocks,  the 
baroclinic  torque  is  nearly  zero;  the  flow  is  nearly  isentropic. 

To  bring  out  the  effects  of  vortex  stretching,  we  consider  the  transport  equation  for  the 
turbulent  enstrophy: 

I 

2~^  ~  ~  +  S{ujiu;2  +  u;30  -  SsjUj)  +  ^i^ijk{p~^  Pkp,p),y  (10) 

The  terms  on  the  right-hand  side  of  eq.  (10)  represent  three  different  mechanisms  of  enstrophy 
variation.  The  first  term  generates  enstrophy  through  the  stretching  of  the  perturbation 
vorticity  lines.  Note  that  in  regions  of  positive  dilatation,  the  generation  of  enstrophy  is 
weakened,  while  it  is  strenghened  when  the  flow  is  compressed.  It  has  been  shown  by 
Betchov  [10]  that  in  incompressible  isotropic  flow,  the  spatial  average  of  this  term  and  (ui,i)^ 
have  the  opposite  sign.  The  velocity  derivative  skewness  is  defined  as 


Sk  = 


(11) 


From  experimental  measurements,  5A:  <  0,  from  which  one  deduces  that  (at  least  for  isotropic 
turbulence)  the  average  of  the  product  of  the  eigenvalues  of  Sij  is  negative  [10].  The  second 
term  in  the  right-hand  side  of  eq.  (10)  represents  the  generation  of  enstrophy  by  interaction 
of  the  mean  shear  with  the  perturbation  flow  field.  Note  the  extra  effect  of  dilatation.  The 
third  term  is  the  baroclinic  torque  combined  with  viscous  effects.  In  general,  if  the  flow  is 
nearly  isentropic,  and  the  Reynolds  number  is  sufficiently  high,  this  term  is  negligible.  We 
do  not  consider  it  further. 


Invoking  a  Helmholtz  decomposition  of  the  velocity  field,  one  can  separate  out  the  effects 
of  compressibility  from  those  solely  due  to  the  intrinsic  solenoidal  turbulent  component  of 
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0 

'I'lie  vortex  stretching  teriii  has  b(‘en  explicitly  rewritten  in  terms  of  and  Sjy  The  first 
coui[-)oneiit  is  now  the  ineoinpn'ssible  vortex  stretcliiag.  The  two  sets  of  numbers  below 
Ivi-  (12)  refer  to  two  simulations.  I'lie  topmost  row  of  numbers  correspond  to  tlie  previous 
run  discussed  (i.e.  S  —  2U),  while  the  lower  set  results  from  an  incompressible  simulation 
with  identical  parameters  exce[)t  that  prms  =  0  and  :\Jt  =  0.  The  total  rate  of  enstrophy 
increase  is  +2,198  for  the  compressible  case,  and  +3,547  for  the  incompressible  run.  Both 
the  prcKluction  terms  and  the  dissipation  terms  are  larger  when  the  flow  is  incompressible. 
Note  that  compressibility  has  little  direct  effect;  all  the  terms  that  involve  the  dilatation 
remain  rather  small  (less  than  1%  of  the  total  amount  of  enstrophy  generation).  However, 
the  compressible  terms  im|)licitly  decreases  the  enstrophy  growth  rate  by  altering  the  balance 
between  vortex  stretching  and  viscous  destruction.  We  also  found  (numerically)  that 
and  are  identically  zero.  I'liis  was  later  proven  rigorously  for  general  homogeneous 

Hows  [24]. 


4  Incompressible  Homogeneous  Shear  Flow 

In  the  past  several  years,  there  have  been  a  series  of  new  findings  related  to  the  small 
scale  structure  of  incompressible  homogeneous  turbulence.  These  findings  are  the  result  of 
investigations  into  the  origin  of  intermittency  in  turbulence.  It  was  found  [25,  16,  27]  that 
there  are  localized  regions  of  very  intense  vorticity  (more  than  4  times  the  mean  value), 
tubular  in  shape,  with  a  length  to  diameter  ratio  of  about  6  [28].  These  results  have  been 
generally  obtained  from  numerical  simulations  with  R\  >  100  which  is  a  microscale  Reynolds 
number  far  in  excess  of  H.\  ^  40  in  our  simulations.  This  might  explain  why  direct  evidence 
of  these  tubes  has  not  been  found  in  our  numerical  data.  However,  the  vorticity  does  align 
itself  with  the  intermediate  principal  direction  of  S,j  when  the  dissipation  is  high.  It  was 
also  discovered  that  the  rate  of  strain  ellipsoid  has  a  very  unique  shape,  with  eigenvalues 
in  the  ratio  (-4:1:3),  in  strongly  dissipative  regions  [13].  'llie  oliserved  tubular  structure 
of  vorticity  is  counter-intuitive,  since  two  positive  eigenvalues  imply  that  vorticity  should 
be  stretched  in  two  extensional  directions  thus  leading  to  sheet-like  patterns.  .+  theory  to 
explain  this  has  been  advanced  by  Majda  [18].  Let  Ai  <  <  A^  be  the  three  eigenvalues 

of  According  to  Majda’s  scenario,  the  vorticity  is  initially  stretched  in  the  direction  of 
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Aa,  the  largest  positive  eigenvalue.  The  evolution  equation  for  neglecting  the  effects  of 
pressure  and  viscosity  then  implies  that  the  rate  of  strain  ellipsoid  is  stretched  in  the  direction 
perpendicular  to  A3.  Thus  the  eigenvalues  Ai  and  A2  increase  in  time  at  the  expense  of  A3. 
At  some  point  the  intermediate  and  expansive  eigenvalues  interchange  their  roles  while  the 
vorticity  retains  its  previous  orientation,  which  is  now  along  the  intermediate  eigenvector. 
While  this  explanation  is  certainly  attractive,  results  by  She  et  al.  [15]  indicates  that  while 
A2  increases  as  predicted  by  Majda,  it  never  exceeds  A3.  An  explanation  must  be  sought 
elsewhere.  A  simplified  model  of  the  Euler  equations  explains  this  alignment;  howev'er,  the 
physical  interpretation  of  the  results  is  unclear. 

The  pressure  gradient  also  shows  preferred  alignment  with  the  principal  directions  of 
Results  of  Ashurst  et  al.  [13]  indicate  that  Vp  is  aligned  with  the  most  compressive  piinci[jal 
direction,  when  the  pressure  gradient  is  high. 

A  better  understanding  of  the  small  scale  dynamics  (and  its  relation  to  coherent  struc¬ 
tures)  will  ultimately  lead  to  better  turbulent  closures,  and  perhaps,  if  we  are  lucky,  to  a 
description  of  turbulence  with  far  fewer  degrees  of  freedom  than  the  standard  scal¬ 

ing.  Any  hope  of  successfully  modeling  compressible  turbulence  thus  hinges,  in  part,  on  an 
understanding  of  how  compressibility  affects  some  of  results  found  in  incojiipressible  flows. 
Are  the  alignment  properties  of  u  and  Vp  still  present  in  compressible  homogeneous  flows? 
Are  these  properties  statistically  the  same  in  regions  of  expansion  (0  >  0)  and  compres¬ 
sion  (0  <  0)?  If  the  flow  properties  are  insensitive  to  compressibility,  there  is  hope  that 
the  methodologies  developed  for  incompressible  flows  will  extend  nicely  to  the  compressible 
cases  with  only  minimal  modification. 


5  Compressible  homogeneous  shear  flow 


In  this  section,  we  summarize  some  of  known  results  of  compressible  homogeneous  shear  flow 
turbulence  that  are  pertinent  to  this  paper.  The  only  extensive  direct  numerical  simulations 
of  compressible  homogeneous  shear  flow  to  date  are  those  of  .Sarkar  et  al.  [1]  and  Blaisdell  et 
al.  [22].  They  cover  different  parameter  regimes,  and  obtain  results  that  are  generally  consis¬ 
tent  with  each  other.  Sarkar  et  al.  [1]  perform  two  series  of  simulations,  differentiated  solely 
by  the  initial  conditions.  In  the  first  set  of  simulations,  the  initial  data  is  incompressible,  i.e. 
{Prinsh  =  0,  .Vo  =  0  while  the  initial  turbulent  Mach  number  (A/()u  i^  varied.  In  the  second 
class,  (A'/()o  =  0.3  while  {prms)o  and  (Xrms)o  are  varied.  In  all  these  runs,  =  ‘21,  and 
{SI\ft^)o  =  7.2  where  I\  =  f'he  turbulent  kinetic  energy  and  is  the  solenoidal 

dissipation.  The  major  result  is  that  increased  levels  of  compressibility,  whether  measured 
by  an  increase  in  Mt,  or  by  an  increase  in  prms  (or  x)  decrea.ses  the  growth  rate  of  kinetic 
energy. 

Sarkar  et  al.  [1,  2]  compute  the  probability  distribution  functions  (pdfs)  of  velocity, 
vorticity  and  dilatation.  As  expected,  the  velocities  (related  to  the  larger  scales)  display  a 
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Gaussian  distribution  (characterized  by  a  skewness  Sk  ^  0,  and  a  flatness  «  3.  The  single 
exception  to  this  was  for  the  transverse  velocity  with  a  slightly  higher  flatness  of  3.3.  Results 
also  indicate  that  the  dilatation  held  is  more  intermittent  than  the  vorticity  held,  and  that 
regions  of  compression  are  more  likely  to  occur  than  regions  of  expansion.  This  is  true  even 
in  the  absence  of  shocks.  No  evidence  of  shocks  were  found  in  these  simulations  because  Rx 
is  too  low.  Blaisdell  [22]  finds  similar  results,  including  at  least  one  flow  with  shocklets. 


6  Methodology  of  determining  statistics 


We  consider  data  from  DNS  on  128^  grids.  A  Helmholtz  decomposition  is  applied  to  the 
velocity  field  and  the  solenoidal  and  irrotational  rate  of  strain  tensors  are  then  computed. 
To  simplify  the  subsequent  analysis  of  the  irrotational  rate  of  strain  tensor,  it  is  made 
traceless  by  subtraction  of  0<5,j/3.  The  derivatives  of  the  velocity  are  calculated  with  a 
order  compact  scheme  rather  than  with  spectral  collocation  to  reduce  internal  code  storage 
requirements.  A  description  of  the  algorithm  can  be  found  in  [29].  Given  the  components  of 
the  rate  of  strain  tensor,  we  compute  its  (real)  eigenvalues  (A,,  i  =  1,2,3)  which  are  ordered 
from  smallest  to  largest:^ 

Ai  <  A2  <  A3.  (13) 

Along  with  the  eigenvalues,  the  three  principal  directions  are  also  computed  and  normalized 
to  unit  length.  The  eigenvector  is  labeled  f,-  with  a  superscript  C  ox  1  depending  on 
the  rate  of  strain  tensor  from  which  it  is  derived.  When  studying  the  orientation  of  these 
unit  vectors,  we  are  only  concerned  with  direction,  not  orientation.  Therefore,  we  take  the 
absolute  value  of  the  direction  cosines  before  calculating  their  distribution. 

In  what  follows,  we  consider  two  types  of  diagnostics.  The  first  type  is  a  simple  one¬ 
dimensional  pdf,  sometimes  weighted  by  a  second  scalar.  The  second  type  of  diagnostic  is 
a  simplified  form  of  conditional  averaging.  Given  two  variables  U\  and  we  define  U\{U2) 
as  a  function  of  IA2-  The  notation  Ui{U2)  should  be  interpreted  as  the  average  of  U\  over  the 
regions  where  U2  lies  within  a  small  band  Ai/2-  For  all  plots,  the  range  of  the  independent 
variable  is  forced  to  lie  within  ±5  rms  from  its  mean  value.  This  range  is  then  divided  into 
25  equal  bins.  These  then  determine  the  discrete  pdf.  All  pdfs  are  normalized  so  that  their 
integral  over  the  independent  variable  is  unity. 

To  reduce  computer  memory  requirements  when  performing  the  diagnostics,  we  only  use 
sample  points  from  the  lower  left-hand  96^  points  in  the  domain.  This  provides  more  than 
800,000  sample  points.  In  the  plots  that  follow,  a  lack  of  sample  points  is  characterized  by 
oscillations  in  the  distributions.  To  clearly  show  the  lack  of  sample  points  (when  it  occurs), 
these  distributions  are  not  smoothed. 

^Note  that  many  authors  order  the  eigenvalues  from  largest  to  smallest. 
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7  Results 


We  consider  results  from  one  simulation  with  initial  data  S  =  20,  M^q  =  0.2,  i/q  —  1/150, 
R\o  =  13.5  and  incompressible  initial  conditions.  Most  of  the  plots  are  shown  at  5/  =  9 
when  Rx  =  23.4  and  AR  —  0.27. 

In  this  section,  we  study  several  properties  of  the  rate  of  strain  tensor  to  help  identify 
characteristics  of  the  flow  which  are  solely  the  result  of  compressibility.  Because  the  flow  is 
anisotropic,  it  is  useful  to  consider  characteristics  of  the  flow  both  with  respect  to  the  fixed 
laboratory  coordinate  system,  and  with  respect  to  a  coordinate  system  attached  to,  and 
rotating  with  an  individual  fluid  element.  Global  orientation  properties  of  are  established 
by  studying  how  the  principal  directions  of  Sij  are  distributed  in  the  laboratory  coordinate 
system.  (In  isotropic  turbulence,  this  diagnostic  is  of  course  devoid  of  meaning.)  Intrinsic 
small  scale  properties,  which  are  hopefully  more  independent  of  the  large  scale  conditions, 
are  characterized  by  considering  the  relative  orientation  of  vectors  such  as  vorticity  and 
scalar  gradients  with  respect  to  the  principal  directions  of  Sij. 


7.1  Correlations 

Based  on  DNS  numerical  databases,  Blaisdell  [22]  showed  that  in  high  shear  homogeneous 
flows,  not  only  are  0^  and  uncorrelaied  which  implies  that 

=  (14) 

but  they  are  also  statistically  independent  which  implies  that 

/(0V)  =  /(e2)iV)  (15) 

where  f{x)  is  the  pdf  of  the  variable  x.  We  confirm  this  result  by  plotting  the  pdf  of  vorticity 
conditioned  on  0  (Figure  1).  It  is  obvious  that  the  distribution  of  vorticity  is  approximately 
independent  of  whether  the  fluid  is  compressing  (0  <  0)  or  expanding  (0  >  0).  This 
implies  that  the  average  of  within  regions  whose  dilatation  is  in  a  specified  range  is 
independent  of  that  range.  Furthermore,  high  vorticity  and  high  dilatation  regions  are 
distributed  independently  of  each  other. 

One  must  be  careful  w'hen  extrapolating  small  correlations  to  zero  correlations.  Unless 
the  correlation  is  exactly  zero,  Eq.  (14)  is  not  strictly  true.  When  two  variables  IF  and  Q  are 
found  to  be  uncorrelated,  it  is  useful  to  consider  the  fraction 

c{rs)  =  (16) 

This  fraction  should  be  very  close  to  unity  for  uncorrelated  variables.  We  find  that  C(0^,»v’^)  = 
0.99,  so  that  Eq.  (14)  is  valid.  Two  other  pairs  of  variables  are  found  to  be  uncorrelated  in 

this  sense.  They  are  (u;^,5'fj  .S^)  and  {Sfj 
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Tlie  eustrophy  convlalcs  ratluM'  strongly  vvilli  'I'able  1  lists  vaiious  pairs  of  vari¬ 

ables,  and  their  correlations.  'The  dilatation  correlated  with  any  variable  constructed  from 
the  soleiioidal  velocity  produces  correlation  coeflicienls  less  than  10“^. 


dabh'  1:  Correlation  factor.s  bet  .oen  several  pairs  of  How  variables 

d  he  most  interesting  feature  of 'Table  1,  is  the  strong  correlation  of  with  -S’/jS/j,  while 
the  streamwise  and  transverse  vorticity  correlate  very  ]>uorly  with  S'ljSjy  On  the  other  haiid, 
w'j-  and  Uy  correlate  stia)ngly  with  each  other,  which  is  consistent  with  a  prefered  v  u’ticity 
orientation  in  the  .r  —  y  pHane. 

7.2  Intrinsic  structure  of  5,; 

.•\s  mentioned  earlier,  in  regions  of  high  c.issijnition,  the  rau  ol  strain  tensor  has  a  preferri'd 
shajie  when  the  How  is  homogeneous  and  incompressible. 

The  soleiioidal  coni])i.;nent  of  retains  its  preferred  sliaije  in  the  presence  of  cc>ii;j>ress- 
ibility.  This  is  shown  in  f  igure  2.  The  pdf  of  A.^/A^  [jeaks  at  a  ratio  of  A.^/A^  =  —0.2. 
When  the  pelf  i-  weighted  with  res|)ect  to  >he  soleiioidal  dissipation  the  peak  of  the 

distribution  is  sharjx'r,  and  has  shifted  to  a  ratio  closer  to  A.^/A^  =  —0.25,  which  corresponds 
to  the  eigenvalues  Aj  in  the  ratio  ( — 1  :  1  :  3).  This  result  matches  the  Hndings  of  .Ashurst  et 
al.  [13].  d  his  eigenvalue  ratio  is  almo'^l  constant  across  the  range  ol  mean  shear  considered 
in  the  1).\S.  .A  case  that  cor; csiioiids  to  isotropic  turbulence  is  presented  in  Figure  4. 

O 

We  now  ask  whether  .V  also  has  some  prelerred  shapi*  in  regions  ol  high  irrotational 

O  O 

strain,  measured  by  the  am|)litude  of  .S' C  ,9  C.  'Phe  results  are  shown  in  Figure  3  where 
t  he  (anit  last  between  weighted  and  unweighted  ilisti  ibiit  iuus  is  much  mure  protiouticed.  '1  he- 
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maximum  value  of  the  weighted  pdf  is  more  than  twice  the  maximum  of  the  unweighted  dis¬ 
tribution,  thus  indicating  that  regions  of  strong  compressible  dissipation  impose  a  structure 
different  than  in  the  rest  of  the  flow.  In  Figure  3,  the  pdf  peaks  for  an  eigenvalue  ratio  of 
(-2.2  :  1  :  1.2). 

We  have  not  been  able  to  establish  a  definite  preferred  ratio  of  A^/Af,  independent  of 
the  mean  shear  rate,  contrary  to  what  is  found  for  A.^/Af.  In  Figure  4,  the  pdfs  of  A.2/A[ 
corresponding  to  a  compressible,  shockless,  isotropic  simulation  (S=0)  clearly  show  that  its 
behavior  with  respect  to  weighting  by  S-jS-j  is  identical  to  that  found  at  5  =  20  (Figure 

0  O 

2).  However,  weighting  of  A^/Af  with  respect  to  Sfj  Sfj  (Figure  5)  does  not  significantly 
affect  its  pdf  when  the  flow  is  isotropic.  Thus,  the  mean  shear,  the  structure  of  the  pressure 
field,  and/or  the  properties  of  the  viscous  forces  must  be  the  major  driving  force  behind  the 

O  ^ 

preferred  structure  of  Sfj.  Further  study  is  clearly  reciuired  before  a  choice  is  made.  .\ote 
that  for  isotropic  turbulence,  exact  solutions  to  approximations  to  Euler's  equations  have 
established  that  the  preferred  shape  of  Sjj  is  due  to  either  the  viscous  stress  terms,  or  the 
properties  of  the  pressure  Hessian  (p.q). 

O 

To  help  understand  the  evolution  of  Sfj,  we  consider  time  histories  of  the  unweighted 
pdfs  of  both  A2/A[  and  A^/Af.  In  Figures  6-7,  we  plot  the  pdfs  of  A^/Af  and  A^/Af  at 
St  =  3,5,9,11,13,21.  After  an  initial  transient  {St  <  9),  the  pdf  of  Aj/Af  reaches  a  quasi 
steady-state  with  a  peak  at  an  eigenvalue  ratio  close  to  -0.25.  The  steady-state  reached 
by  Aj/Af  is  not  as  definitive  as  that  of  Aj/Af  (Figure  7)  judged  on  the  amount  of  scatter 
between  different  times.  Furthermore,  A^/Af  reaches  a  quasi  steady-state  at  an  earlier  time 
{St  =  5)  than  does  Aj/Af.  This  slight  difference  of  timescales  will  resurface  when  we  consider 
the  orientation  of  Sij  with  respect  to  the  laboratory  coordinate  system. 


7.3  Extrinsic  Orientation  of  5,y 

By  definition,  in  isotropic  turbulence,  the  principal  directions  of  the  rate  of  strain  tensor 
are  randomly  oriented  in  the  laboratory  coordinate  system.  However,  homogeneous  shear 
turbulence  is  anisotropic,  and  therefore  we  expect  the  rate  of  strain  ellipsoid  to  have  a 
preferential  orientation.  This  is  obvious,  since  the  mean  shear  itself  imposes  a  constant 
strain  on  the  flow,  oriented  at  45°  to  the  x  and  y  coordinate  directions. 

Consider  the  angle  Oij  between  eigenvector  f,  and  the  unit  coordinate  axis  vector  x,  from 
which  the  pdfs  of  |  cos  Oj^l  and  |  cos  are  computed.  The  cosine  of  Oij  is  chosen  rather  than 
Oij  so  that  the  probability  density  function  of  a  Gaussian  distribution  curve  is  fiat.  Because 
the  independent  variable  is  the  cosine  of  an  angle,  one  must  be  careful  when  interpreting 
the  magnitude  of  the  pdf  maxima.  Strong  peaks  of  the  |cosf^|  pdf  do  not  necessarily  imply 
a  strong  peak  of  the  pdf  of  0.  The  angle  at  which  this  peak  occurs  must  also  be  taken 
into  account.  For  example,  in  Figure  8b  the  pdf  reaches  values  of  1.25  in  regions  where 
I  cos  01  ~  0.7  {0  w  45°),  while  in  F'igure  10b,  it  reaches  values  of  4  for  |  cos0|  s;  1.0  {0  ~  0°). 
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Since  d&  =  —d{cos  0)j  sin(6^),  a  given  d{cos  6)  corresponds  to  a  much  larger  dO  a.i  6  =  0°  than 
at  =  45°. 


Figures  8a-c  show  the  pdfs  of  |cos^/j|  at  St  =  9.  In  this  and  the  next  figure,  each 
subplot  corresponds  to  one  eigenfunction.  It  is  evident  that  the  directions  of  maximum 
extensional  and  maximum  compressive  strain  have  qualitatively  similar  distributions.  These 
two  eigenvectors  have  alignments  which  are  most  often  at  angles  near  45°  to  the  x  and  y 
axes.  However,  the  45°  orientation  with  respect  to  the  streamwise  direction  is  strongest.  This 
suggests  that  the  mean  strain  plays  the  role  of  a  polari2er  which  forces  these  two  directions 
to  align  with  itself.  The  pdf  of  j  cos  02jl  shows  a  tendency  to  align  at  30°  to  the  x  axis.  When 
the  pdfs  of  ^  are  weighted  with  respect  to  S-jS-j,  f/  and  fF  show  an  increased  tendency  to 
align  at  45°  to  the  x  and  y  axes  as  seen  in  Figure  9.  On  the  other  hand,  the  weighted  pdf 
of  I  cos02i|  indicates  that  f2^  has  equal  probability  of  lying  at  53°  to  each  of  the  3  axes.  Note 
that  this  alignment  property  of  ij  is  much  weaker  than  that  of  f/  and  as  determined  by 
the  ratio  of  maximum  to  minimum  probability. 

The  pdfs  of  |cos61F|  are  markedly  different  than  those  of  |cos^/^|  as  seen  from  a  com¬ 
parison  of  Figures  8a-c  and  Figures  lOa-c.  Once  again,  the  most  compressive  and  the  most 
expansive  principal  directions  have  almost  identical  pdfs.  The  strongest  trend  is  that  and 
align  with  the  y  axis. 

The  intermediate  irrotational  eigenvector  behaves  differently.  This  eigenvector  orients 

itself  primarily  along  the  x  and  z  axes.  It  is  also  evident,  that  fp  has  a  strong  propensity 

to  remain  normal  to  the  y  axis.  The  analysis  has  been  repeated  using  data  from  several 

turbulent  simulations  at  Mt  =  0.2, 0.4,  shear  levels  of  S’  =  10,20,30,  and  most  of  the 

general  conclusions  remain  unchanged.  The  only  noticeable  effect  of  weighting  the  pdfs  of 
0  0 

I  cos6J,^|  with  Sfj  Sfj  is  to  sharpen  the  maxima  (Figure  11).  The  qualitative  features  of  the 

distributions  are  not  affected  by  the  weighting.  In  summary,  the  most  extensional  and  most 

compressive  principal  directions  of  S-j  tend  to  be  at  45°  to  the  coordinate  directions,  while 

0 

the  principal  directions  of  S'^  tend  to  align  with  the  coordinate  directions. 

Why  are  the  properties  of  Sfj  so  different  from  those  of  S-j?  Although  the  dynamics  of 

0 

S’/j  are  somewhat  understood,  this  is  not  the  case  for  Sfj.  When  the  flow  is  incompressible, 
the  vortex  stretching  term  in  the  enstrophy  transport  equation  is  of  the  form  iOiS-jUj.  The 
vorticity  stretching  is  initially  in  the  direction  of  the  maximum  mean  strain,  i.e.  at  45°  to 
the  streamwise  direction.  This  in  turn  initially  aligns  the  principal  directions  of  S-j  so  that 
fg  is  aligned  with  w.  However,  the  change  in  uj  induces  a  change  in  the  vorticity,  whose 
interaction  with  S-j  serves  to  reorient  S-j  relative  to  w.  This  process  is  described  by  the 
analytic  formulations  of  Vieillefosse  [19],  Majda  [18],  and  Cantwell  [17],  which  predict  that 

O 

the  vorticity  will  eventually  align  itself  with  The  evolution  of  Sf^,  on  the  other  hand,  is 
driven  by  one  of  two  major  effects: 

1 .  The  structure  of  5,^  acts  as  a  source  term  or  as  spatially  dependent  coefficients  in  linear 


12 


equations  for  Sfj,  or 

O 

2.  The  mean  shear  is  primarily  responsible  for  the  structure  of  Sfj. 

To  investigate  which  of  the  two  scenarios  is  the  most  likely  candidate,  we  plot  the  dis¬ 
tributions  of  [cos^ijl  at  St  =  3,5,9,11,13,21  (Figures  12a-c  to  Figures  15a-c)  Each  plot 
shows  the  pdf  (at  various  times)  of  one  of  the  9  directional  cosines.  From  Figure  12,  we 
find  that  after  an  initial  transient,  the  pdfs  of  |cos0(,|  tend  towards  a  stationary  shape.  As 
expected,  these  distributions  are  mirrored  by  those  of  |  cos^3,  |.  On  the  other  hand,  it  is  not 
clear  that  all  the  components  of  |cos02il  reach  an  equilibrium  shape  during  the  simulation 
(see  Figure  13).  For  example,  the  orientation  of  with  respect  to  x  was  still  changing  when 

the  numerical  simulation  ended.  The  next  two  figures  show  the  equivalent  distributions  for 

0 

ff  (Figures  14a-c  and  Figures  15a-c).  The  existence  of  a  stationary  structure  of  Sf  is  much 
more  definite,  and  occurs  earlier  than  the  equilibration  of  f/.  Therefore,  it  is  likely  that  the 
mean  shear  rather  than  Sjj  determines  the  compressible  rate  of  strain  tensor  to  zeroth  order. 


7.4  Orientation  of  the  vorticity  vector 

As  expected  from  previous  investigations  [13, 15]  the  vorticity  orients  itself  along  in  regions 
of  high  solenoidal  dissipation.  If  the  angle  between  u  and  f/  is  denoted  by  this  trend 
is  clearly  seen  by  conditioning  cos^0,-  on  S-jSfj  (Figure  16).  In  regions  of  high  dissipation, 
u  aligns  itself  strongly  in  the  direction  of  f2,  and  normal  to  both  fi,  and  fa.  This  result 
is  consistent  with  results  of  Ashurst  et  al.  [13]  and  She  et  al.  [15].  From  Figure  8b,  fa  is 
oriented  at  approximately  30°  —  45°  at  St  =  9,  which  is  consistent  with  the  results  of  Rogers 
and  Moin  [31].  They  find  that  as  the  shear  rate  increases,  the  angle  that  the  vorticity  forms 
with  the  streamwise  direction  decreases  from  45°  to  approximately  30°. 

The  compressible  dissipation  has  no  influence  on  the  relative  orientation  of  uj  and  f/. 
Indeed,  the  pdfs  of  |cos^^/‘^|  are  flat  when  conditioned  by  0,  as  seen  in  Figure  17.  This  is 
consistent  with  a  lack  of  correlation  between  the  dilatation  and  the  vorticity.  The  vorticity 
and  f^  show  a  weak  tendency  to  align  with  one  another  (Figure  18).  Although  not  shown 
here,  we  have  found  that  the  average  alignment  of  the  vorticity  vector  with  f/  is  essentially 

independent  of  either  0  or  Sf  Sf. 


7.5  Orientation  of  the  pressure  gradient  vector 

Also  of  interest  are  the  properties  of  the  pressure  gradient.  Referring  back  to  Eq.  (6)  for  the 
transport  of  0  and  its  higher  moments,  it  is  clear  that  the  properties  of  the  pressure  gradient 
and  pressure  Hessian  must  be  understood  before  there  is  any  hope  of  correctly  modeling 
the  subtle  relationships  between  dilatation  and  vorticity.  As  shown  in  Erlebacher  et  al.  [21], 
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tlie  total  pressure  (in  isotropic  turbulent  flows)  can  be  thought  of  as  a  combination  of  an 
elliptical  (incompressible)  component  which  together  with  satisfies  the  incompressible 
Navier-Stokes  equations,  while  the  difference  between  the  pressure  and  i.e.  the  com¬ 
pressible  pressure  ,  is  representative  of  compressible  effects.  An  alternate  decomposition 
of  the  pressure  field  is  considered  in  Sarkar  [32],  in  which  the  incompressible  pressure  re¬ 
tains  all  the  effects  which  invoh'e  the  mean  density  gradients.  This  generalization  of  the 
pressure  decomposition  is  found  useful  for  the  analysis  of  inhomogeneous  turbulent  flows. 
Both  decompositions  have  similar  characteristics  regarding  the  properties  of  the  pressure 
components.  The  incompressible  pressure  is  basically  elliptic,  i.e.  information  in  one  part  of 
the  flow  is  immediately  transmitted  throughout  the  flowfield.  In  contrast,  the  compressible 
pressure,  p^  is  hyperbolic,  i.e.  the  [)ressure  and  dilatational  fields  are  responsible  for  waves 
traveling  with  the  speed  of  sound  in  the  flow. 

The  presence  of  both  compressible  and  incompressible  character  in  the  pressure  field  can 
be  at  least  partially  extracted  from  the  data  by  considering  the  orientation  of  the  pressure 
gradient  with  respect  to  f/  and  .  In  Figure  19,  we  plot  the  pdf  of  the  pressure  gradient 
alignment  with  respect  to  f/.  W'e  observe  that  there  is  a  preferred  orientation  of  Vp  at  about 
-10°  to  both  f/  and  f/.  In  fact,  the  pdfs  of  |cos<l[‘^j  and  |cosf^3^|  are  almost  identical,  sug¬ 
gesting  that  with  regard  to  the  pressure  gradient,  the  most  extensive  and  most  compressive 
principal  directions  of  the  solenoidal  rate  of  strain  tensor  play  similar  roles.  .Ashurst  et  al.  [13] 
find  that  in  incompressible  homogeneous  shear  flow  turbulence,  Vp  tends  aligns  itself  in  the 
dizection  of  f/.  This  contrasts  with  our  findings  which  indicate  that  Vp  is  inclined  similarly 
with  respect  to  f/  and  In  I'igure  19b  we  have  plotted  the  pdf  of  Vp^  The  greatest 
difference  with  Figure  19a  if  the  decreased  sensitivity  of  the  pdf  of  Vp^  to  orientation. 

In  Figure  20,  we  show  the  variation  of  the  average  of  |cosf^/^P  as  a  function  of  0.  The 
orientation  of  the  pressure  gradient  with  respect  to  f/  shows  only  a  slight  dependence  on 
dilatation. 

Figure  21  shows  the  pdfs  of  [  cos  dj"  ^|.  The  pressure  gradient  is  aligned  with  or  f;p  with 
equal  [zrobability.  In  addition,  there  is  a  strong  influence  of  the  compressible  dissipation  on 

the  alignment  of  with  ff  (Figures  22-23).  Botli  5^  and  0  are  used  to  condition  the 

distributions  of  cos^fl;  .  In  regions  of  high  {Sfj  Sfj),  Vp  is  strongly  aligned  with  fi,  while 
it  has  quite  different  orientations  in  regions  of  strong  expansion  and  strong  compression. 
When  0  >  0,  the  alignment  is  towards  f^,  while  it  is  in  the  direction  of  when  0  <  0. 
,\s  expected,  the  solenoidal  component  of  pressure,  p\  lesponds  differently  to  the  dilatation 
(Figure  24).  I'he  orientation  of  Vp^  with  respect  to  f/  shows  much  weaker  dependence 
on  with  0  than  does  Vp.  fhe  gradient  of  p'  has  different  mean  orientation  in  regions  of 
compression  and  expansion.  When  0  <  0,  Vp^  is  perpendicular  to  ,  while  it  is  somewhat 
aligned  with  in  regions  of  flow  expansion.  For  completeness,  we  also  show  in  Figure  25 
the  orientation  properties  of  V p‘  with  respect  to  f/  when  conditioned  on  S-jSly  Although 
the  correlation  is  not  strong,  in  regions  of  strong  solenoidal  dissipation,  Vp^  is  pulled  slightly 
towards  the  intermediate  eigenvector.  On  the  other  hand,  Vp  shows  a  prefered  alignment 
with  ff . 
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From  Figures  22-23,  one  also  infers  that  Sfj  Sfj  and  0  are  correlated  (see  Table  1).  This  is 

corroborated  by  Figure  26  where  the  average  magnitude  of  dilatation  decreases  monotonically 
0  0 

as  Sfj  Sfj  is  increased.  This  type  of  plot  can  easily  be  misinterpreted  however,  as  evidenced 

O  O 

by  Figure  27  which  shows  that  the  average  of  Sf  Sf  is  maximum  in  regions  of  high  10|.  The 
bias  of  0  in  regions  of  high  compressible  dissipation  results  from  the  skewness  of  pdf(0). 


8  CONCLUSIONS 


Compressible  turbulence  is  a  difficult  subject  due  to  the  number  of  parameters  which  control 
the  flow.  The  purpose  of  this  work  was  to  determine  what  are  the  leading  order  effects  that 
compressibility  has  on  the  small  spatial  scales  in  homogeneous  compressible  shear  flow  turbu¬ 
lence.  To  this  end,  the  logical  candidates  to  study  are  the  rate  of  strain  tensor,  the  vorticity, 
the  dilatation  and  the  pressure  gradient.  In  addition,  the  already  successful  Helmholt:/  de¬ 
composition  was  put  to  good  use  to  further  separate  out  solenoidal  from  irrotational  behavior 
in  the  flow. 

The  data  analysis  was  done  by  considering  the  invariants  of  S,j.  Of  crucial  importance 
were  its  eigenvalues.  The  orientation  of  its  eigenvectors,  both  with  respect  to  a  coordinate 
system  fixed  in  the  laboratory  and  with  respect  to  w  and  Vp,  were  studied.  The  effect  of 
dilatation  and  of  solenoidal  dissipation  were  considered  using  conditional  sampling. 

We  find  that  the  qualitative  statistical  properties  of  the  vorticity  and  solenoidal  rate  of 
strain  tensor  are  affected  very  little  by  compressibility  effects.  This  is  due  in  part  to  the  lack 
of  correlation  between  the  dilatation  and  the  vorticity  fields.  The  only  major  compressibility 
effect  is  on  the  orientation  of  the  pressure  gradient  In  regions  where  the  flow  is  compressed, 
Vp  aligns  with  ff ,  while  it  is  oriented  along  in  expansion  regions.  This  effect  is  also  seen 
in  compressible  isotropic  turbulence,  but  is  much  weaker. 

O 

In  shear  flow,  the  eigenvectors  of  S-j  and  Sf  have  a  preferred  orientation  with  respect 
to  a  frame  of  reference  fixed  in  the  laboratory.  While  the  principal  directions  of  5^  tend  to 
orient  themselves  at  angles  of  30°  -  50°  to  each  of  the  three  coordinate  directions  (with  the 

O 

strongest  alignments  along  the  x  axis),  those  of  Sf  are  primarily  aligned  along  the  axes  of 
laboratory  coordinate  system.  Both  tensors  have  an  orientation  which  tends  to  a  stationary 
distribution  for  large  time.  However,  the  irrotational  rate  of  strain  tensor  ‘"equilibrates” 
on  a  faster  time  scale,  which  suggests  that  it  is  the  mean  shear  that  is  responsible  for  its 
properties,  not  the  fluctuating  rate  of  strain  S'j. 

While  there  are  an  infinite  number  of  ways  to  combine  variables,  an  equal  number  of 
ways  to  display  them,  not  to  mention  interpret  what  is  displayed,  we  feel  that  our  initial 
attempt  has  succeeded  in  opening  a  veritable  Pandora’s  box  of  new  questions  to  be  answered. 
For  example,  are  regions  of  strong  vorticity  and  strong  dilatation  correlated  at  a  distance? 
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Perhaps,  regions  of  strong  vorticity,  together  with  the  mean  shear,  induce  the  irrotational 
velocity  field  to  acquire  a  particular  structure.  What  are  the  roles  played  by  and  p^ 
in  the  dynamics  of  these  structures?  Can  the  structure  of  the  irrotational  rate  of  strain 
tensor  be  explained  by  linear  theory,  or  are  nonlinear  effects  the  primary  driving  force? 
These  are  but  some  of  the  questions  which  come  to  mind.  It  will  be  interesting  to  witness 
how  the  information  presented  herein  will  eventually  get  blended  into  a  coherent  theoretical 
framework. 
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Figure  2:  Pdf  of  Aj/Af  with  and  without  weighting.  S  =  20,  u  =  1/150, 
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Figure  3:  Pdf  of  A^/Af  with  and  without  weighting.  5’  =  20,  i/  =  1/150 
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Figure  5:  Pdf  of  Aj /Af  with  and  witl 


pdf  (x^Af) 


Figure  6:  Pdf  of  weighting)  at  5f  =  3, 5, 9, 11, 13, 21.  5  =  0,  i/  =  1/150,  St  =  9. 
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Figure  7:  Pdf  of  A^/Af  (no  weighting)  at  St  =  3,5,9,11,13,21.  5  =  0,  i/  =  1/150,  St  =  9. 
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P'igure  11;  Pdf  of  jcos j  =  (1,2,3)  weighted  by  Sfj  S^.  Mt  =  0.2,  5  =  20,  =  1/150, 

St  =  9.  (a)  i  =  1,  (b)  i  =  2,  (c)  i  =  3. 
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Figure  12;  Pdf  of  |cosf^(,|.  .\/;  =  U.2,  S  =  20,  u  =  1/150,  St  =  3,5,9,11,13,21.  (a)  i 
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Figure  13:  Pdf  of  |  cos02^|-  =  0.2,  5  =  20,  t/ =  1/150,  Sf  =  3, 5, 9, 11, 13, 21.  (a)? 
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Figure  14:  Pdf  of  ]  cos  ^f,|.  Mt  =  0.2,  5  =  20,  i/ =  1/150,  =  3,5,9, 11, 13,21.  (a) 
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Pdf  of  |cos<^^;|.  .\/,  =  0.2,  S  =  20,  u  =  1/150,  St  =  3,5,9,11,13,21.  (a)  i 
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Figure  16:  Variation  of  leos^t^/-*]  with  //,  conditioned  on  S-^Sj 
=  1/150,  St  =  9. 
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Figure  19a:  Pdf  of  [cos^rn.  Mt  =  0.2,  5  =  20,  =  1/150,  bt 


Figure  19b:  Pdf  of  |  cos  0/^*’  |.  A/,  =  0.2,  5  =  20,  v  =  1/150,  St  =  9. 


Figure  26:  Variation  of  0  conditioned  on  Sfj  Sfj.  Mt  =  0.2,  S  =  20,  u  =  1/150,  St  =  9. 


Figure  27:  Variation  of  Sfj  Sf  conditioned  on  0.  Mt  =  0.2,  5  =  20,  i/  =  1/150,  =  9. 


37 


REPORT  DOCUMENTATION  PAGE 

form  Approved 

OMB  A/o  0704'0f8a 

Public  feoo77i''-  Dyrae'’  rnis  oi'eaicri  of  ''’C"ra;  ot'  'S  *o  .i.e'd-?**  *  f'our  oe^  'e^porse  including  tne  lime  to^  reviewing  mstruaions.  searcning  e«  stmg  aata  sources 

gathering  joC  maintaining  tne  aata  needed,  and  -  ocTDieting  an©  revi^^win©  rre  ruenion  ot  i.nformation  Send  comments  regarding  this  Durden  estimate  or  arty  otrter  asDec*  ot  this 
collection  ot  intormation.  nciuOmg  suggest 'Ons  ‘  '  reducing  tris  Duroen  t:  Aiasnington  Headquarters  Services.  Direaorate  for  mtot^mation  Ooerations  and  Reoons.  12  iS  vetferson 
Davis  Highway,  Suite  ’2C4  ..ningion  uA  2  2202-1304.  and  \c,  tne  O’tice  ot  Manaqemerst  and  Budget  Paoerwor*-  Reduction  Project  {0704-0188)  Washington.  DC  20503 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

Anri] 1  Q9? _ r.nnf rar rnr  Rpnnrr _ 

4.  title  and  subtitle 

STATISTICAL  ANALYSIS  OF  THE  RATE  OF  STRAIN  TENSOR  IN 
COMPRESSIBLE  HOMOGENEOUS  TURBULENCE 

5.  FUNDING  NUMBERS 

C  NASl-18605 

WU  505-90-52-01 

6.  author(S) 

G.  Erlebacher 

S.  Sarkar 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Institute  for  Computer  Applications  in  Science 
and  Engineering 

Mail  Stop  132C,  NASA  Langley  Research  Center 

Hampton,  VA  23665-5225 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

ICASE  Report  No.  92-18 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

National  Aeronautics  and  Space  Administration 

Langley  Research  Center 

Hampton,  VA  23665-5225 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  N"VRER 

NASA  CR- 189640 

ICASE  Report  No.  92-18 

t1.  SUPPLEMENTARY  NOTES 

Langley  Technical  Monitor:  Michael  F.  Card  Submitted  to  Physics  of  Fluids 

Final  Report  A 

12a.  DISTRIBUTION/ AVAILABILITY  STATEMENT 

Unclassified  -  Unlimited 

Subject  Category  34 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

Recent  analysis  of  direct  numerical  simulations  of  compressible  homogeneous  shear 
flow  turbulence  has  unraveled  some  of  the  energy  transfer  mechanisms  responsible  for 
the  decrease  of  kinetic  energy  production  when  the  flow  becomes  more  compressible. 

In  this  complementary  study,  we  focus  our  attention  on  the  rate  of  strain  tensor. 

A  Helmholtz  decomposition  of  the  velocity  field  leads  us  to  consider  a  solenoldal  and 
an  Irrotatlonal  rate  of  strain  tensor.  Their  eigenvalue  distributions,  eigenvector 
orientations,  and  the  relative  alignment  between  the  eigenvectors  and  the  vorticity 
and  pressure  gradient  vectors  are  examined  with  the  use  of  probability  density  func¬ 
tions.  The  Irrotatlonal  rate  of  strain  tensor  is  found  to  have  a  preferred  structure 
in  regions  of  strong  dilatation.  This  structure  depends  on  the  mean  shear,  and  is 
quite  different  from  that  of  solenoldal  rate  of  strain  tensor.  Compressibility 
strongly  affects  the  orientation  properties  of  the  pressure  gradient  vector. 

14.  SUBJECT  TERMS 

compressible;  turbulence;  rate  of  strain;  alignment;  statistics 

IS.  NUMBER  OF  PAGES 

39 

16.  PRICE  CODE 

A03 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

NSN  7540-0 1-280- 5500  Standard  Form  298  (Rev  2-89) 


Pfe\<rtD^  by  AN^’  Std  Z)9-’9 
299-102 


NASA  l.anRlrv,  l«92 


